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The assessed phase diagram of the boron-carbon system contains a single non- 

bJ[^ stoichiometric boron-carbide phase of rhombohedral symmetry with a broad, ther- 

| modynamically improbable, low temperature composition range. We combine first 

principles total energy calculations with phenomenological thermodynamic modeling 

to propose a revised low temperature phase diagram that contains two boron-carbide 
■ i— i ■ 

^ | phases of differing symmetries and compositions. One structure has composition B4C 

and consists of BnC icosahedra and C-B-C chains, with the placement of carbon 
on the icosahedron breaking rhombohedral symmetry. This phase is destabilized 
above 600K by the configurational entropy of alternate carbon substitutions. The 
\ other structure, of ideal composition B13C2, has a broad composition range at high 

rH ■ 

temperature, with rhombohedral symmetry throughout, as observed experimentally. 
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I. INTRODUCTION 



The phase diagram of the boron-carbon system is highly controversial and has been 
frequently revised. In addition to the pure elements, anywhere from one to eight different 
compound phases are claimed l|, |2j, with the sole point of agreement being the existence of 
a nonstoichiometric boron-rich rhombohedral phase known colloquially as "boron carbide" . 
Boron carbide is a hard material that is useful in making armor. Owing to the neutron 
absorbance of boron, it also has uses as a shielding, control and shutdown material in 
nuclear power plants. In the assessed phase diagram reflecting current consensus, reproduced 
in Fig. HI a single phase, labeled as B 4 C and exhibiting rhombohedral symmetry, covers a 
composition range from a low carbon content of 9% to a high of 19.2%. Note that the 20% 
carbon content implied by the substance name B4C is never achieved. 

The broad composition range indicates a substitutional solid solution. The temperature 
independence of the phase boundaries makes the assessed phase diagram thermodynamically 
improbable Specifically, if the phase field were to extend to absolute zero (T=0K), 
the positive entropy associated with substitutional disorder would imply a violation of the 
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FIG. 1. Assessed phase digram of the boron carbide system 
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third law of thermodynamics [5[. Several examples of such apparent violations are known, 
and generally are resolved by the onset of new phase behavior at low temperatures. One 
famous example is the dispute over the phase diagram of plutonium-gallium jfj], where the 
American phase diagram extending the composition range of the 5-phase to low temperature 
was eventually rejected in favor of the Russian version in which S is stable only at elevated 
temperature. 

Very likely, boron carbide is out of thermodynamic equilibrium at all but the very highest 
temperatures. In this case, conventional experimentation cannot easily determine the equi- 
librium phase diagram, and theoretical approaches are needed to resolve the fundamental 
character of the stable state(s) at low temperature. Previous first principles calculations 
utilizing electronic density functional theory jl-lO] indicate the existence of two boron-rich 
phases, one of ideal stoichiometry B13C2 and rhombohedral symmetry, the other of ideal 
stoichiometry B 4 C=B 12 C3 and monoclinic symmetry. In the present paper, we introduce 
a simplified thermodynamic model, inspired by first principles calculations, that elucidates 



the probable evolution of the equili 
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ibrium boron carbide phase field at low temperatures. 
jj-[l3[ claim that the B 4 C phase has a 15-atom unit cell 
with rhombohedral symmetry. These papers agree that the boron carbide structure can be 
viewed in terms of principal structural elements: boron icosahedra linked in a fashion similar 
to the a-boron structure, and a linear chain of 3 atoms lying along the 3-fold axis of the unit 
cell. The atomic positions of the icosahedra belong to two site classes: polar sites linking 
the icosahedra to each other, and equatorial sites linking the icosahedron to the chains. 

The distribution of carbon and boron atoms on these structural elements is uncertain. An 
idealized B 13 C2 structure with rhombohedral symmetry occupies the twelve icosahedral sites 
with boron, while the three-atom chain takes the pattern C-B-C (i.e. the chain center atom 
is boron and the terminal sites are carbon). However, first principles calculations HQ 
find that the most stable, enthalpy minimizing, structure is B12C3, with one carbon atom 
replacing boron on a polar site of the icosahedron, breaking the rhombohedral symmetry. 
Experimentally |15|], the phase appears to contain a mixture of characteristic motifs: B12 
and B n C icosahedra; C-B-C, B-B-C, B-V-B and C-V-C chains (V = vacancy). 

The enthalpy minimizing B12C3 structure has monoclinic symmetry, in disagreement 
with the B 4 C phase's observed rhombohedral symmetry. We propose that the phase known 
as "B 4 C" should be renamed "B 13 C2", or simply "rhombohedral", as we shall do for the 
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remainder of this paper. We also claim that there exists a second phase, which is the 
true B 4 C phase, that we term "monoclinic" . In the monoclinic phase, all of the carbons 
lie on equivalent polar sites of the icosahedra, making a well-ordered structure. Although 
there is an entropic term corresponding to the choice of this site, it is non-extensive and 
thus vanishes in the thermodynamic limit. In the absence of substitutional disorder the 
monoclinic phase is a line compound. In previous work [loj], we proposed a phase transition 
in which the monoclinic phase transitions into the rhombohedral phase through the unlocking 
of a rotational degree of freedom in the placement of the polar carbons. Note that swapping 
of carbon sites is equivalent to a rotation of a BuC icosahedron. Landau theory 16J] predicts 
no such transition occurs in the rhombohedral phase, as the symmetry already matches the 
high T limit. Here, we construct an analytic model to interpret our previous computational 
results, from which we can derive an actual phase diagram. 



II. FREE ENERGY MODEL 

Four solid phases compete for stability: elemental boron ((3- rhombohedral), elemental 
carbon (graphite), and monoclinic and rhombohedral boron carbide. The first three phases 
are modeled as line compounds throughout the entire temperature range, but rhombohedral 
boron carbide will be allowed a carbon composition ranging from B 14 Ci (x = 1/15 = 0.067) 
to B 12 C3 (x = 3/15 = 0.200). The stable T=0K composition for the rhombohedral phase 



is B13C2 (x = 2/15 = 0.133). According to first principles calculations [10] the most stable 
B14C1 structure consists of B12 icosahedra with B-B-C chains. At composition B12C3 the 
most stable structure is monoclinic (i.e., not rhombohedral). Rhombohedral structures at 
this composition correspond to B n C icosahedra with C-B-C chains, but the placement of 
the carbon atom on the polar sites is randomly oriented among different primitive cells, 
unlike the monoclinic structure where all polar carbons are uniquely aligned throughout 
space. Thus our model for the rhombohedral phase allows as structural units B12 and B11C 
icosahedra, and C-B-C and B-B-C chains. Let yc be the fraction of icosahedra containing 
a polar carbon, and ys the fraction of chains containing a terminal B. Note that ys and yc 
are bounded between and 1, and that the carbon fraction 

x= ^(yc-VB + 2). (l) 
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The enthalpies of formation of /3-rhombohedral boron and graphite vanish by definition, 
while we denote the T=0K enthalpies of the stable rhombohedral and monoclinic boron 
carbide phases, respectively, as h° R and h° M . Our first principles calculations yielded h° R = 
—0.087 and h° M = —0.117 eV/atom (see Fig. HI])- Notice that the monoclinic phase is 
more stable, at x — 0.200, than the rhombohedral phase is at x = 0.133. We extend the 
rhombohedral phase entropy beyond its ideal composition by assigning an enthalpic penalty 
> for each excess chain boron and an enthalpic benefit 7 < for each icosahedral polar 
carbon. Thus h(y B ) = h R + (3y B provided yc = and similarly h(y c ) = h° R + jyc provided 
y B = 0, as illustrated in Fig. [Ill In general, y B and yc are both nonzero, so that 

h R (x) = h° R + f3y B + jy c . (2) 

To complete the free energy model we assign an entropy based on random substitution 
disorder. Specifically, an ideal site substitution entropy —k B (ylny + (1 — y) ln(l — y) for 
selecting the fraction y of structural units on which to substitute, together with an intrinsic 
rotational entropy k B In 6 for the orientation of the icosahedral polar carbon, and an intrinsic 
reflection entropy k B \n2 for the choice of terminal chain position, as there are 6 polar sites 
per icosahedron and 2 terminal sites per chain. Thus, 

k 

sr(x) = -^(y B ln2 - y B kiy B - (1 - y B ) ln(l - y B )) 
15 

k B 

+ — (y c ln6 - y c ^yc - (1 -yc)ln(l - yc)) (3) 
15 

combines with the enthalpy eq. (TSJ) yielding the free energy 

g R (x, T) = h R (x) - Ts R (x). (4) 



A similar free energy model was proposed by Emin [17J , although he supplemented the free 
energy with an additional bipolaron density that is no longer considered as relevant. As the 
monoclinic phase is a line compound, no composition dependence is required for its enthalpy, 
and because the rotational degree of freedom is locked, it lacks entropy. 

Fig. [Ill illustrates components of the free energy. Black circles denote the enthalpies h° R 
and h° M , while the black line segments are the convex hull of free energy g(x) at T = OK. Red 
lines and green curves illustrate enthalpy and T = 2000K free energy under the constraint 
that either y B or y c vanishes, taking parameter values /5=0.05 and 7 = —0.024. The blue 
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FIG. 2. Components of free energy model. Details in text. 

curve illustrates the T = 2000K free energy gn(x) (eq. (jlj)) with the constraint relaxed, while 
the dashed pink line is the full convex hull of free energy g(x). 

The free energy gu{x) in eq. (jl]) is to be regarded as a Landau-type free energy, as the 
order parameters y B and yc must still be determined as a function of composition x and 
temperature T. Note the two order parameters are not independent, as the composition 
eq. ([I]) implies (dyc/dy B )\ x = 1. Now, minimizing g R with respect to yc yields 

0=%^ = p + i + W(-l nl2 + ln( * B / g )) (5) 
dye x 15 V. \1 -y B l -ycj J 

Substituting for the parameter y B and rearranging yields a quadratic equation for yc, 

y 2 c (l - k) + (2 - 15x(l - K))y c + k(1 - 15x) = (6) 



where we define 



k(T) = i2e- 15 ^/ kBT 



(7) 



7 



as a measure of the extent to which polar boron atoms can swap positions with chain terminal 
carbons, which constitutes the fundamental excitation of the ideal rhombohedral structure. 

A. T ^ OK limit 

We require that the enthalpy parameters + 7 > to ensure that the ideal rhombohedral 
structure with ys = yc — minimize the free energy at x = 2/15 in the limit T = 0. Hence, 

lim k(T) = 0, (8) 

T->0 

reducing eq. (jHJ) to the equation y c + (2 — 15x)yc = with two solutions: either yc = 
(and then ys = 2 — 15x), or else yc = 15a; — 2 (and then ys — 0). That is, either boron 
substitutes on chain terminal carbon sites, in which case yc = and x < 2/15, or else 
carbon substitutes on icosahedral polar boron sites in which case ys = and x > 2/15. As 
a function of composition we have 

j/ B = (2-15x)0(2/15-x), j/c = (15z - 2)0(z - 2/15) (9) 

where is the Heaviside step function. 

Although exact only at T = OK, the essential singularity in k makes this an excellent 
approximation over a wide range of temperature, up to 1000K and sometimes even higher 
in the examples discussed later. The approximation represented by eq. ([9]) separates the 
free energy into two branches, an x < 2/15 branch where only ys terms contribute, and an 
x > 2/15 piece where only yc terms contribute, as illustrated in Fig. [Til In each case, one 
branch of this piecewise-analytic free energy competes with a line compound, so we next 
derive a general equation for the phase boundary in this scenario. 

1. Substitutional Disorder Coexisting with a Line Compound 

To find the phase boundaries of rhombohedral boron carbide we must locate the coex- 
istence of our substitutional^ disordered phase with the competing phases. Depending on 
composition and temperature the coexisting phase might be /3-rhombohedral boron, mon- 
oclinic boron carbide, or graphite. In every case the coexisting phase is treated as a line 
compound whose free energy is simply its enthalpy. In this section we solve the coexistence 
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equations generally, then apply this solution to specific phase boundaries in the following 
sections. 

Consider the free energy model for substituting on a fraction y of structural units, each 
with intrinsic multiplicity Q and enthalpic penalty 5, 

l rp 

g(y,T) = h° R + 5y-^(y\nn-y\ny-(l-y) ln(l - y). (10) 

Such a free energy represents one of the two branches of our rhombohedral free energy 
(eq. (|4])) at low temperature. Let this phase coexist with an ordered line compound of free 
energy g* = h* and "composition" y*. We now wish to find the composition, y', of the 
disordered phase that coexists with the ordered phase at y*. 

Coexistence is determined by a double tangent condition, at y' and y*. Specifically, there 
exists a straight line f(y) — fo + fiy that is tangent to g* at y = y* and to g(y,T) at 
y = y' . For the line compound, tangency at y* is the simple condition f(y*) = h*. For the 
disordered phase, tangency requires both that f(y') = g(y',T) and that f' = fi= g'(y',T), 
where 

g\y,T) = 5-^(\nn-\n T ^-y (11) 

Solving, we find 

n v* e iB(ft*-hi-v)/* fl T = ( y >y*(i - y ') l ~y\ (12) 

In the case where y' ~ 0, this equation simplifies to 

y' « fl e i6(fc*-h&-W**kBT (13) 

while for y' m 1, we have 

y' « 1 - (n e iHh*-h° R -6y*)/k B T^ ~t-y* (14) 

2. Boron-rich phase boundary 

We apply this general solution to the specific case of rhombohedral boron carbide coex- 
isting with /3-rhombohedral boron. In the notation of the preceding section, h* = 0. There 
are two cases to consider depending on whether the phase boundary lies to the left or to the 
right of x = 2/15. 

Case 1. If the boundary lies at x < 2/15, so that yc = 0, then we identify y = ys, 
5 = (3 and Q = 2, as the disorder corresponds to substitution of boron onto the terminal 
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chain carbon sites. Also, y* = 2 corresponds to the composition x = 0. The requirement 
that /3-rhombohedral boron be stable at x = against the boron carbide phase at ys = 2 
implies a constraint that 2/3 > — h° R . The phase boundary occurs in the limit of small y, and 
from eq. ( fl3|) we have 

y' B = 2e -15(^+2/3)/2*BT_ (15) 

Case 2. If the boundary lies at x > 2/15, so that y R = 0, then we identify y = yc, 
5 = 7 and Q = 6, as the disorder corresponds to substitution of carbon onto the icosahedral 
polar sites. Also, y* = —2 corresponds to the composition x = 0. The requirement that 
free energy be convex at x = 2/15 implies a constraint that 27 > h° R . The phase boundary 
occurs in the limit of small y, and from eq. ( fl3|) we have 

y' c = 6e +15(/ l f A-27)/2fc S T_ 

To determine if case 1 or 2 occurs, consider the ratio 

y' c /y' B = 3e-^- h °^ T . (17) 

If 7 — (5 — h° R is positive, then the ratio vanishes at low temperature and the requirement of 
convexity places us in case 1. If instead it is negative, then the ratio diverges and we have 
case 2. 



3. Carbon rich phase boundaries 

Next we apply our general solution to the specific case of rhombohedral boron carbide 
coexisting with monoclinic boron carbide. Now, in the general notation, h* = h° M . As in 
the case of coexistence with elemental boron, two specific cases are possible. However, the 
fact that carbon substitution is known to be energetically favorable (i.e. 7 < 0) implies 
that only the case of ys — and y = yc is relevant. Thus we identify 5 = 7 and Q = 6. 
Noting that the condition for low temperature stability of monoclinic B 4 C against disordered 
rhombohedral atx = 3/15is/ije + 7> h° M , the phase boundary occurs in the limit of small 
y', and from eq. ( fl3l) we have 

y' c = 6 e -^(h° R -h M +7)/k B T (lg) 

as the phase boundary in coexistence with the monoclinic phase. 
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However, above a certain temperature T the rhombohedral phase coexists with graphite. 
We assume T lies in the low temperature limit, and find the general form for the phase 
boundary due to coexistence between graphite and the rhombohedral phase. In our general 
notation, y* = 13, y = yc, £ = 7 and Q = 6. In contrast to the preceding cases, the phase 
boundary occurs in the limit of y m 1, and from eq. (114jl we have 

,15(^ + 13/3), , s 

y' c = l- 6-->exp( y ^ ,J> ) (19) 

as the phase boundary in coexistence with graphite. The two carbon-rich boundaries cross 
at a certain temperature T that can be determined by setting the values of y' c equal in 
eqs. ( 1T8]) and ( TT91 . For realistic parameters, the crossing occurs at yc ~ 1. Then from 
eq. ( ITS]) we find 

T w 15(h° R -h° M + j)/k B In 6. (20) 



B. Analytic T=oo Limit 

We now examine the high temperature limit for the free energy model. For all possi- 
ble values of (3 and 7 we have lim^oo k(T) = 12 giving the quadratic equation —\\y c + 
(2 + 165x)yc + 12(1 — 15x) = 0, only one of whose roots is physical: yc = ^(2 + 165x — 



V532 — 7260x + 27225a; 2 ). Notice that only the multiplicities enter into this equation. Al- 
though this equation is more complicated than the one obtained for the low temperature 
limit, it is defined across the composition range from x = 1/15 to x — 3/15. The phase 
boundary on the boron rich side, x' boron , is given by 

dgn . h/3-boron 9 r( x borons 

T) 



Arp ' boron C\ rpl 

adj u x boron 



(21) 



and the phase boundary on the carbon rich side x' carbon is given by 

dg R h 

graphite 9 R\ x carbons 

T) 



rjiy ' carbon 1 rp? 

ax 1 x carbon 



(22) 



The high temperature limits are dominated by the entropic terms, yielding nonlinear implicit 
equations for x' . Using a numeric equation solver we find x' boron = 0.1095 and x' carbon = 
0.1515. However these values are reached only at extreme high temperatures, while our 
model is intended only for use below the melting temperature. 
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Parameter 


Value (eV/atom) 




-0.0869 


h° M 


-0.1167 


a 


0.1031 


7 


-0.0244 



TABLE I. Parameters for the rhombohedral phase obtained from first principles calculations. 

III. REALISTIC PARAMETER VALUES 

Our simple model depends on just four parameters, h° R , h° M , (3 and 7. Of these, the values 
o f h R. 

M are easily determined from first principles calculations with simple idealized models. 
Estimated values of and 7 may be obtained by inserting a single B or C substitutional 
defect into a hexagonal supercell of the ideal rhombohedral B 13 C 2 structure. A general 
description of the computational method is in Ref. [h]]. Resulting values are listed in 
Table IHIi These parameters obey the constraints discussed in previous sections. Because 
7 — — h° R = —0.0406 is negative, we are in case 2 as discussed in Section III A 21 where 
x ' Boron > 2/15 at low temperature, although it eventually goes to x' Boron < 2/15 at very high 
temperatures. 



There is some question whether our structural model is complete in the B-rich limit, as 
alternate structure models contain additional sites in the chain region, some only partially 



occupied |13|, |l_8|, Il9j . To take into account a possible influence of these additional sites, 



we investigate the effect of reducing the value of /3, in order to model the effect of lower 
boron-rich enthalpies. Note that increased multiplicity Q would also enter the free energy 
linearly in y R , though with an added factor of temperature T. 

Results for a selection of values of /3 are shown in Fig. IHII As expected, the monoclinic 
phase is destabilized above a temperature To ~ 600K, at which it decomposes into a co- 
existence of carbon-rich rhombohedral phase together with graphite. This value of To is 
surprisingly consistent with the location of the heat capacity peak previously reported that 
was obtained from a completely different method [10]. The maximum carbon content of the 
rhombohedral phase is bounded below 20%, owing to the logarithmic singularity in s(y) at 
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FIG. 3. Predicted phase boundaries for our model using calculated parameters as listed in Table ITTT1 
and a selection of values of f3. Line compounds are /3-boron at x = 0, graphite (not shown) at x = 1 
and monoclinic boron carbide at x = 3/15. Carbon-rich phase boundaries of rhombohedral boron 
carbide are shown as solid lines, while boron-rich phase boundaries are shown as dashed lines. 

y c = 1 creating an infinite slope in g(x,T) (to weak to be visible in Fig. [TTJ) . Note that the 



Gibbs phase rule ^| implies that rhombohedral and monoclinic boron carbide must have 
differing compositions while in coexistence with graphite, hence x' c = 0.2000 is forbidden in 
principle. 

At temperature 2400K (around the melting point) the maximum carbon content depends 
on the value of /3, and ranges from 19.9% down to 19.7% in our model for the given range 
of f3 considered, whereas the experimentally assessed limit is 19.2%. The phase boundary 
in coexistence with boron depends strongly on the value of 0, as can be seen in Fig. IIHI 
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The experimentally assessed limit is 9% carbon, but over the range of (3 values studied here, 
the limit ranges from 10-16% carbon, strictly above the assessed value. This is a further 
indication of inadequacy of our model in the boron-rich limit. 



IV. CONCLUSION 



In this paper we propose a simple free energy model for boron carbide. We predict the 
existence of two low temperature compounds, resolving the mystery of the assessed low 
temperature composition range. One phase, whose ideal composition is -B13C2, has rhombo- 
hedral symmetry throughout its wide high temperature composition range, consistent with 
experimental observations. The other phase is a line compound of exact stoichiometry B 4 C 
characterized by an array of parallel BnC icosahedra whose symmetry is monoclinic. De- 
spite its favorable low enthalpy, this phase is predicted to be stable only below T = 600K, 
a temperature so low that it might not be possible to form in thermodynamic equilibrium. 
Various limiting behaviors of the model phase boundaries are derived analytically, including 
the key fact that x' Carbon < 0.200 in the rhombohedral phase. A quantitative discrepancy 
in the boron-rich limit reveals the need for further enhancement of our model. Effects to 
consider include the additional interstitial sites, as well as electronic and vibrational {20 1 
entropy. 
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